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This part of the theory directly continues the derivation of the kinetic 
equation presented in the preprint 990958 of this archive. All definitions are 
the same and no special remarks are given. 

1 Chapman-Enskog procedure 

At first we recall the standard version of the Chapman-Enskog procedure 
and then we shall present the generalization to our situation and fulfill cal- 
culations. 

Consider equation 

f = -(^ + B)/> (1) 

where operators A and B have the following properties 

• Operator A depends on variable x, operator B depends on x and y. 
Both operators don't depend on time t. 

• Operator B is supposed to be small in comparison with A. Later we 
shall rescale time as to have ||y4|| ~ 1, ||i?|| ^ 1. The value ||i?|| will 
be the small parameter of the theory. 

• The eigenvalues and eigenfunctions of B are unknown. The eigenfunc- 
tions Ai and eigenvalues Oj of operator A are supposed to be known. 
They have the following properties 

Qq = 



ai> 6 > i ^0 



with some positive small delta. The eigenfunctions are supposed to be 
normalized. 

One can introduce symbol 0(| |) which means that the given expression 
has the order of or less. 

Initial conditions for (|l|) are supposed to be known. Then one can directly 
describe the evolution during the first periods of time. Here one can neglect 
B. One has only to solve equation 

dP 

Due to the known eigenfunctions of A this solution is quite obvious 

oo 

P = exp{-amt)PmAm (3) 

m=0 

where 

Pm = {P\t=o, A^) (4) 

will be called the mode amplitudes {Am will be modes) and the symbol ( , ) 
denotes the scalar product. 

The final of relaxation according to is obvious 

P^PoAo (5) 
The characteristic time of relaxation is 



trel = l/Ol (6) 

where oi is the smallest eigenvalue of A (except the zero value). 

Now A{pqAq) = and operator A exhausts its power. One has to consider 
the action of 

^ The physical reason to consider the further stages of evolution is that in the nucleation 
problem Aq — I and the final state of relaxation doesn't provide ant benefit into the value of 
the nucleation rate. It lies in the evident correspondence with the formula J ^ exp{Fc) / 
of the one-dimensional theory where stands the big parameter Av. 
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The deviation of P from PqAq can be caused only by the action of operator 
B. So, in some sense it is smalL Then we shall seek the solution in the 
following form 

P = poAo + A (7) 

For the correction term A in accordance with the perturbation theory one 
can consider the representation 

oo 

A = E5^^^ (8) 

i 

where every term S^^^ has the order 

5^'^ Oi\\B\\') 

Now we are going to consider the time dependence of the amplitudes Pm- 
One can note that the action of A leads to the vanishing of pm for m ^ 0. 
The unique source to fill the given mode is to take into account the action 
of operator B which is small but still can provide the transition from one 
mode to another. But B doesn't depend on time. So, the intensities of 
transitions between modes don't depend on time also. So, the unique source 
of time dependence is the time dependence of the zero mode po- Later this 
dependence will be spread by B to all other modes. 

As the result we state that the time dependence of an arbitrary mode 
amplitude is going through the time dependence of po: 

Pmif) = Pm{.Poit)) (9) 
dprn ^ dprndpo , 

dt ~ dpo dt ^ ^ 

Certainly the time dependence of po isn't fixed and can not be obtained 
immediately. One has to recalculate it at the every step of approximation^. 
The value S^^^ depends on time t through po{t) in the i-th. approximation. 
At every newQ 5^*^ one has new p^ and new Pq\ So, one has a new time 
dependence. 

^It is also possible to leave this dependence as some formal parameter and to calculate 
it when the "final" approximation (i.e. the approximation with the necessary precision) 
is obtained. 

•^The lower index denotes modes. 
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The substitution of expansion 



p = po^ + E^^'^ (11) 



i=l 

into (|1]) gives 

+ E ^ = -{A + B)ip,Ao + E 5«) (12) 

Note that the i-th term has order i, but the value with the order i will have 
benefits from all terms with j > i. 

The r.h.s. and l.h.s. of the last equation are some functions and we sup- 
pose that they can be decomposed into series over the small parameter ||-B||. 
We shall denote the i-th term which has the order 1 1 -B 1 1 * by the superscript 
< i >. Then one can write for the r.h.s. 

oo 

{{A + B){poAo + E5«))<*> = A<5« + B5^'~'^ (13) 

i=l 

The consideration of the l.h.s. is more complicated. At first we have to 
see the expression for dpo/dt. Having projected equation (|l|) on Aq one can 
see that 

^ = -(Ao,5P) (14) 
Hence, one can come to the estimate 

^ = 0(||i3||) 

One can rewrite the last relation using the decomposition (^ 

^ = -(Ao, 5(^0^0 + E5«)) (15) 

Now it is evident that dpo / dt contains all powers of the small parameter and 
it has the order 0(||i?||). Now it is worth to determine the j-th term in 
decomposition over parameter It quite easy to do for dpo/dt. Namely, 

(^)<^> = -(Ao,i?5(-^)) (16) 
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Then one can fulfill the same transformations for dd^'^ydt. One present 
the following estimate 

So, one can reduce the infinite sum to the finite sum 

(E^r>-(E^r> (18) 

The next transformations are evident 

(E^)- = (E|^t)-= (19) 
-E^(^o,s<5(^-^-^)) 

Now we can substitute all results into the initial equation and get 

- (Ao, B5^'-'^)Ao - E ^(^0, B5^'-^-'^) = (20) 

-A6^^ - B5^'-^'^ 

As far as the eigenfunctions of A are known one can use the last relation 
to get 5^^. 



2 Relaxation 

Consider now our situation. Certainly, one can choose operator Di as the 
main one. 

An operator 

^ ^ , d 
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has the known set of eigenfunctions which are the Hermite polynomials Hj 
(j = 0, 1, 2, 3, 4, ). The scalar product is defined as 



1 /■°° 

JtX J-oo 



which leads to 

The eigenvalues of Dq are 



X- = -2i i = 0,1,2, 



One has to mention two important relations 



d_ 

The first one gives the way to construct the set of Hermite polynomials. 



Both these equations allow to introduce in |IU| the representation of modes. 
We shall call the operators in the r.h.s. of two previous equations as the 
transition operators. The first one is the mode increase operator. 

The second one is the mode decrease operator. Then one can consider 
the operators of multiplication on fi and the differentiation over /x as the 
superposition of transition operators. Any additional operator appears from 
the non-Fokker-Planck behavior and from higher terms in the free energy. 
Non-Fokker-Planck terms give high derivatives. Free energy can be well ap- 
proximated in the nearcritical region by polynomial and leads to additional 
multiplications on /i. So, the additional term can be regarded as the super- 
position of the mode increase operators and mode decrease operators. This 
procedure resembles the formalism of the secondary quantization. Below we 
shall use another approach. 

It is clear that Di can be presented as 



1=1 
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with coefficients 



= - E ^^S^.H'+I m = 2,3,4, 



4" =E + 

i' ^ H'^j j 

Then the eigenfunctions of Di will be the Hermite polynomials^ and the 
eigenvalues are given by 

A,=i:4"(--'— ^' — 



1=1 



As the result it can be seen that one can take Di as the main operator in 
the Chapman-Enskog procedure. The approximate form of kinetic equation 
will be the following one 

But already in the investigation of the stationary distribution establishing one 
will see that the formal priority of Di in comparison with D2 isn't sufficient 
to ensure the relaxation to the stationary state. 

Operator D2 will be small in comparison with Di in terms of some small 
parameter (not in formal sense described above) only when tj far all j are 
small parameters. Then the main terms in operator D2 will be the terms 
with the smallest sum of indexes I + m. Then the required condition of the 
smallness of D2 formulated in terms of the relaxation times will be 



*In the main order S2j goes to 1. 



7 



This parameter has been required^ to be small in [^]. Then the solution at 
the relaxation stage can be written as the seria of the relaxation modes 

oo 
1=0 

where 

are initial amplitudes. This seria ensures the relaxation to the stationary 
state Pq'^Hq-. 

P ^ pJ°)/7o = Prel 

The time of relaxation trei is given by 

trel Di = Ai 

The mentioned restriction prei < 1 (actually according to |^ one has 
to require Prei -C 1 ) practically excludes the nucleation under the strong 
thermal effects. So, the situation considered in p is rather poor. 

One can note that the final state Prei corresponds to the equilibrium 
distribution over yU. So, it is the eigenfunction of the operator Di + D2 with 
the zero eigenvalue 

(D, + D2)Prel = 

This can be easily proven when we return to the form of the finite differences 
and reconstruct Di and D2. 

Certainly it is difficult to determine all eigenfunctions and eigenvalues of 
Di + D2. We needn't all relaxation modes but only the final one. Also we 
have to estimate the time of relaxation. One can prove that the relaxation 
time trel D1+D2 for equation 

dP 

is less or equal to the previous relaxation time trei Di 

trel D1+D2 ^ trel Di 
""with account of the mentioned error. 
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Operator D2 "helps" the relaxation to Prei- The way to prove the last es- 
timate is to consider the blocks along fj, axis and estimate the action of D2 
between blocks. 

As far as D^, have small parameter in comparison with Di then no 
special condition is required. Now we can consider Di + D2 as the main 
operator at the relaxation stage. The smallness of Prei isn't now required. 



3 Correction terms 

When P is close to Prei which is the eigenfunction with the zero eigenvalue 
the operators Di or Di + D2 can not be considered as the main ones. One 
has to use the Chapman- Enskog asymptotic decomposition. But to fulfill 
this decomposition one has to know all eigenfunctions of the main operator. 
So we can not use Di + D2 as the main operator now. 

We shall redefine the main operator after the end of the relaxation stage. 
Now the main operator will be Di. One can see that "correction operator" 
-D2 + -D3 + -D4 has a complicate structure. Operators and have small 
parameter. Operator D2 is small only in the formal sense. So, the natural 
modification of the Chapman-Enskog procedure is to consider two Chapman- 
Enskog procedures. One can include the whole internal Chapman-Enskog 
procedure into every step of the external Chapman-Enskog procedure. The 
initial procedure has the aim to "invert" the action of operator D2 and the 
external procedure will "invert" the action of -(- D4. 

Every new approximation will contain the small parameter 1/Ak ~ k^/^ 
( or at least Ak/k ~ in the order of the number of external approxi- 

mation. This order decrease the order of flow calculated in the initial (zero) 
approximation. But one has to stop the external procedure until the result- 
ing order is greater than = 1 because the continuous description along i/j 
will be violated in this order. One has to note that the evolution along i/j 
is principally different in comparison with evolution along /i. The domain 
along Ui is concentrated only at the integer values. The value of fi is the 
real value. This difference leads to the specific lattice corrections described 



in These corrections have very complicated structure, an account can 
not be fulfilled explicitly in the analytic way. As the result we have to re- 
strict ourselves by the first correction and adopt that the account of the first 
correction in the external procedure is already sufficient. But in the internal 
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procedure we have to calculate the great number of approximations. 

Also one has to note here the requirement to have no operators 5*1^, S@j, 
S^j in the final results. Really, the presence of this operators corresponds to 
the elementary shift taken into account. This lies in contradiction with the 
lattice domain along z/j. 

The matter under discussion is the rate of nucleation and to determine 
this rate and need only the embryos flows averaged over /x. This operation 
corresponds to the projection on Hq ~ const. 

Now we can turn ourselves to establish equations for dpo/dt. The last 
dependence is extremely important because in the Chapman-Enskog proce- 
dure all dependence of approximations on time is going through dependence 
of po on time. To get dpo/dt one can project kinetic equation 

dP 

— = [Di + D2 + D3 + D4\P 
at 

oylHqi.q. fulfill the integration /f^exp(— /x^) dji. The l.h.s. gives dpo/dt. 

Then 

^ = {Ho, [Di + D2 + D3 + D,]P) 
The r.h.s. will lead to more complicate expression. Consider 

{Ho, D,P) = {Ho, D, Y.P^Hi) = {Ho, E ^^P^H^) = AqPo = 

i i 

as far as Aq = 0. Consider 

oo 

{Ho,D2P)^{Ho,D2j2piHi) 

^From explicit form of D2 one can see that the last elementary operator in 
action of D2 will be the mode increase operator — 2jj,). So, if we present 

D2J:ZoP^H^e^ 

QO 00 

D2^PiHi = ^CiHi 

i=0 i=0 

then Co = 0. There is no projection on the zero mode. Then 

00 

{Ho,J2ciHi)^0 

i=l 
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for arbitrary Cj. Then 

{Ho,D2P) = 

Consider {Hq^D^P). Operator can be split into two parts D^^ and 
D^i, defined as 

Dsa = EE - 2^,yLJS^JP{M,^^) 

Then 

(//o,i^3aP)=0 

and the reasons are the same as for D2. The action of D^j, can be presented 
as 

a 00 _m am 

{Ho, Ds,P) - E ^^3.(^0, ^+ E ±i7^P) 



. dfj ' ^ ml dn"^' 



The last operator D4 gives 



As the result one can get 

The values in the r.h.s. of the last equation can be interpreted as the flows 
Jj along Vj axis 

J, = {H,,[L, + WIS,,Y.^^^Q^]P) 

which transfers the last equation in the standard form 

^--V— J- (21) 
dt ~ ^du/' ^^^^ 
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4 Calculations 

Now we can turn to the direct calculations. The distribution P can be 
presented into the following form 

oo 

P = PoHo + Y.PiiPo)Hi 

1=0 
oo 

PliPo) = J2pt'\po) 

m=l 

where the lower index indicates the number of mode and the upper index in 
brackets indicates the number of approximation. As far as we already stated 
we can fulfill in the external procedure only one step and it isn't necessary to 
mark it. The upper index, thus, corresponds to the internal approximations. 

Now we can rewrite expression for Jj in terms of announced decomposi- 
tion. Then 

oo r^rn oo 

J, = (Ho, [L, + w;ss, E fr^iboi^o + j:piipo)m 

m=l 1=0 

Consider various terms of the last expression. Namely, 

{Ho, LjP) = Ljpo 

rn 

Then consider (i^o, Wp^, E^^, ^^[p,Ho+EZoPiiPo)m The termpo^^o 

T™ am 

doesn't leads to any influence b ecause in X^m=i mT 9/*"' there is at least one 
operator of the mode decrease. Then 

oo om oo oo am oo 

(Ho,W;S,j E ±q;^\PoHo + T.Pi(Po)Hi]) = {Ho,W;Ss, E ^qJ^EpMHi) = 

m=l 1=0 m=l "1^ 1=0 

oo oo a oo 

(^0, w;s,, E E ^^PiiPoW - w;ssj E rrPmiPo)2- 

i=l TO=1 m=l 

The resulting expression for Jj will be the following 

oc 

Jj = Ljpo + W;Ssj E r^pM^"" 

m=l 
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In the main order we can omit Ssj. Certainly, the correction term in S^j has 
the order greater than Lj has but there is no zero mode except LjPo. So, we 
have to keep Lj and throw away S^j here. Then 

oo 
m=l 

Now wc can calculate the set of corrections. The initial condition for the 
asymptotic expansion will be the result of the relaxation stage, i.e. Prei- One 
can easy note that 

DlPrel = 
D2Prel = 
DsbPrel — 

because the fist elementary operator is the mode decrease operator. In D^^ 
there are only mode increase operators and then 



°° (—r v 8 
j 1=1 '^A* 

oo (_ \l oo ( \l 

E LJS^, E Si^ (-i)Voi^/ = E L^s,, E ^-jrP^H, 

j 1=1 j 1=1 

The last operator contains the small parameter in the first order. In the main 
order one can throw out. 

Operator has at least the small parameter r in power 2. So, it is small 
in comparison with D^a which has r in the first power. As far as operators 
Di, D2 have no action we have to take into account but we can neglect 
the action of D4. 

As the result we have 
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In the zero approximation 
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The structure of first approximation is rather simple due to the zero 
approximation which is localized only at the zero mode. Now the current 
approximation has all modes and operator D2 will lead to the non-zero result. 
So, the operator D2 will be the main one and one can neglect and D4. In 
D2, thus, one can neglect 

The time derivative can be also neglected. Really, due to ( p^) 

- + + D,]H,) 

in the main order. 

One can easy see that 

{Ho,D2Ho) = Q 
(Ho, DstHo) = 

as far as the first operator in the action of D2 and D^h is the derivative over 
H which gives zero in application to Hq ~ const. 

As far D^a has the small parameter and D4 has the result {dpo/dt)^^^ is 
small. As far as we calculate in the main order of small parameter we can 
neglect {dpo/dtY^\ 

In all approximations from the second one we can use ( PT| ) and see due 
to the smallness of d/dui the smallness of {dpo/dtY^^ i = 2,3, 4, 

We see that the Chapman- Enskog procedure is now reduced to the trivial 
equation (except the zero mode) 

= D^'D2P^'^ 

with the evident initial approximation. In the main order one can take 
away from the operator D2. This form can be obtained from the simple 
analysis without Chapman- Enskog formalism. 

We can present the last equation in the following form 

00 

<?=! 

as the matrix representation of the linear operator Di^D2 in the basis 
(without the zero mode). All can be put to zero. 

One can note the following important features: 1). Any operator acts un- 
til the current moment in the already fulfilled part of the external Chapman- 
Enskog procedure only one time. 2). Every operator has its own specific 
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structure and can not be reproduced by the actions of other operators. Both 
these features allow to forget about Sij, S2j, S^j. If any operator had been 
used more than one time then one would take these operators into account. 

The absence of 5*1^, S2j, S^j is very important in the context of the lattice 
structure of the domain in the {z/j} plane. Only the absence of these operators 
allows to ignore this lattice structure. 



5 Final expressions 

Now we have to estabhsh the expression for Tkq. As far as D2 is the super- 
position of the mode decrease operators and mode increase operators it will 
transfer Hi into superposition of Hi. The action of D^^ on Hj (except zero 
mode) is evident - it is multiplication Hj on \J^. 

How D2 transfers Hq into if^? At first according to the definition of D2 
the mode decrease operators will act. The number I of the mode decrease 
operators will be between Iq and q. Parameter Iq appeared from the evident 
requirement that we have to fall lower than the k-th level. So, Iq = 1 ii 
k > q and Iq = q — k + 1 if q > k. The result of the action of the mode 
decrease operators wiU give the coefficient 2^q\/{q — /)!. Then to get Hk one 
has to apply k — {q — I) mode increase operators which will give coefficient 
(_l)fc-(9-0. As the result 

Also one can consider another representation of T^q. At first we fall from 
the q level to the / level which gives (g^_f)! ff2^~'- then we increase the mode 
from the / level to the k level which gives j^zjy The result will be 

min{q,k}=l q-l , fnti'^~^ 



j ,.„ (9-0!« (*-')! 

Note that F^g doesn't depend on the number of approximation. This will 
lead to some important consequences. The first one is the possibility to write 
equation between p^: 

00 

Pk = J2^kqPq+Pk^ 
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or 

oo 



q=l ^kK- j 

Due to the linearity of equations the decomposition of the initial approx- 
imation will be reproduced in all approximations. Namely, we shall present 
in the following form 

pi" =2:/'!]' 

j 

where 

1 

and the second lower index indicates component in the initial approximation 
(later F^g will mixture different components). 

The linearity results in the possibility of decomposition 



j 



Index j has no correspondence here with the direct decomposition over com- 
ponents. 

For the following recurrent expression 

oo 
g=l 

is valid. 

The total amplitude Pk can be also decomposed as 



Pk = J2Pkj 



where pkj are 



For Pkj one can write 



Pkj^^Pk, 



(i) 
1=1 



.(1) 
kj 

q=l 



Pkj = J2'^kqPqj+Pkj 
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One can also move operators Lj through f = {Tkq}- This gives the seria 
of relations. For the first approximation one can write 

j 

where 

The recurrent relations for af^ will be 



a 



kj — XI ^kqO'qj 
q=l 

In terms of a^*- the value p^'"* can be easily expressed as 



Having introduced 



one can get 



oo 



13 
i=l 



^kj — XI ^kqaqj + a^kj 
q=l 

In terms of akj the value pk can be easily expressed as 

j 

The given decompositions are rather attractive from the first point of 
view, but actually one can not simplify the problem with the help of these 
decompositions. 

The equation for po will the following one 
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and in the main order 



^'^ j '-"^j i 171=1 

where indexes i and j marks components. The last equation has the standard 
form investigated in 0. The part of ^ concerning the solution of the last 
equation in the two dimensional sense is also acceptable as far as it reproduces 
even in details^. 

Now we shall present another method to calculate correction terms in the 
Chapman- Enskog procedure. 

to 



One can easily note that recurrent equations for p\^^ and for will lead 



hm p\^\ —>■ consti A'^. 
lim pf^ consti{'-f, 



max ) 



where 7max is the eigenvalue of f = {Fp^} with the maximal absolute value. 

Then one can say that starting from some number the tails of sums 
HiPf^j^ J2iPt^ resemble the tail of geometric progressions. Then one can 
easily calculate these sums 

E{i) Pi j 

4—m, I max 

^ — "''hm 

and 

oo (jf^lim) 

j—mi I max 

^ — '"'hm 

The first mum terms have to be calculated explicitly. The boundary 
which depends on j, I can be found as the characteristic boundary when 
pj+i^pi approaches some constant independent on i. This procedure gives 
also the value of •jmax- 

Certainly, the most interesting situation is the strong manifestation of 
the thermal effects. Here the thermal effects can not be considered as some 
corrections but radically change the character of the process. 



^except several misprints appeared in Also the initial matrix (with dimension 2) 
isn't diagonal which doesn't produce any difhculties. 
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Really, the main dependence of the nucleation rate is accumulated in 
exponent of the critical embryo free energy. The giant renormalization due 
to the thermal effects means that the temperature of the critical embryo 
really differs from the temperature of the vapor-gas media. The relative 
difference of the temperature expressed in the units of (3 (in estimates we 
can forget about different components) can attain several units. It means 
that the characteristic value of /i attains several units. 

One can easily note that to reproduce the real solution (i.e. to attain the 
values of /x in several units) we have to take into account a great number of 
modes. Really, the Hermite polynomials Hn are the polynomials of power 
n. The characteristic region of localization of the function if„(/i) exp(— /x^) 
is [~ —y/n,r^ -y/n]. So, to describe the situation correctly we have to take 
into account at least /Zq modes. This quantity equals to the dimension of the 
matrix F which becomes also great. This produces numerical difficulties. 

Consider the matrix F. Due to factorials in the denominators the limit 
of the elements with a big indexes is going to zero. It takes place under 
the arbitrary Tj. The r.h.s. of the matrix equation, i.e. the known vector 
= {pf^^ also has vanishing elements when i oo. The structure of the 
matrix F is the following: the upper triangle matrix has the elements with 
big values due to the big coefficient appeared from the action of the mode 
decrease. The lower triangle matrix is "smaller" than the upper one. Thus, 
the following method will be rather effectively one can split matrix F into 
the upper triangle matrix (with big elements) and the lower triangle matrix 
(with moderate elements). Then at every step of the iteration procedure the 
upper triangle matrix will be inverted (it is easy to do). The zero mode (in 
fact it is the first line in the matrix equation) will be calculated explicitly. 

One can easily note the the iteration solutions of the matrix equation 
return us to the initial formulation of the problem at the level of the recurrent 
relations. So, the presented approach to calculate the maximal eigenvalue 
gets now a solid ground. It is also more effective because it allows to estimate 
the necessary number of modes taken into account directly (by the smallness 
of the higher mode at every step). 

^With the special account of the Chapman-Enskog specific features 
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